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ABSTRACT 



Context. Star formation (SF) takes place in unusual places such as way out in the intergalactic medium out of material expelled from parent 
galaxies. 

i-Q Aims. Whether SF proceeds in this specific environment in a similar way than in galactic disks is the question we wish to answer. Particularly, 
we address the reliability of ultraviolet, Ho- and mid-infrared as tracers of SF in the intergalactic medium. 
q ■ Methods. We have carried out a multiwavelength analysis of the interacting system NGC 5291, which is remarkable for its extended HI ring 
$H ' hosting numerous intergalactic HII regions. We combined new ultraviolet (GALEX) observations with archival Ha, 8 pm (Spitzer Space 
, Telescope) and HI (VLA B-array) images of the system. 

■ Results. We have found that the morphology of the star forming regions, as traced by the ultraviolet, Ha, and 8 pm emission is similar. The 
8.0 pm infrared emission, normalised to emission from dust at 4.5 pm, which is known to be dominated by PAH bands, is comparable to 

• i-H , the integrated emission of dwarf galaxies of the same metallicity and to the emission of individual HII regions in spirals. The 8.0 pm in the 

■ intergalactic environment is therefore an estimator of the star formation rate which is as reliable in that extreme environment as it is for spirals. 
There is a clear excess of ultraviolet emission compared to individual HII regions in spirals, i.e. the [8.0] / [NUV] and [Ha] / [NUV] SFR 
ratios are on average low although there are some large variations from one region to another, which cannot be explained by variations of the 
metallicity or the dust extinction along the HI structure. Comparing the observed SFR with a model of the evolution of [Ha] / [NUV] with 
time favours young, quasi-instantaneous though already fading starbursts. The total star formation rate measured in the intergalactic medium 
(which accounts for 80% of the total) surrounding NGC 5291 is up to 1.3 Mq yr~' - a value typical for spirals - assuming the standard SFR 
calibrations are valid. The SFR drops by a factor of 2 to 4 in case the star formation is indeed quasi-instantaneous. 
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1. Introduction 

It has long been known that "vast and often very irregular 
swarms of stars and other matter exist in the spaces between the 
conventional spiral, elliptical, and irregular galaxies" (Zwicky 
195 1). Those stars have been stripped off from their host galax- 
ies during interactions, the intergalactic medium being too dif- 
fuse and hot to allow in situ star formation. However a fraction 
of the stellar population could nevertheless have been formed 
locally out of the gas which has been pulled out from peripheral 
regions of gas rich galaxies. This has been shown to occur in 



simulations of tidal encounters (Hibbard & Mihos 1995; Due 
et al. 2004) as well as in the case of collisional rings (Appleton 
& Struck-Marcell 1996; Horellou & Combes 2001). This gas 
which has the metallicity of the parent galaxy (roughly solar 
metallicity in the case of normal star forming spirals 1 , lower 
for dwarf galaxies, see Mendes de Oliveira et al. (2006) for 
instance) can then be reprocessed to form new stars and even 
new dwarf-like galaxies if the expelled gas reservoir was mas- 
sive enough. How star formation (SF) proceeds in collisional 
debris between galaxies, in an environment which is isolated 
from the direct influence of the galactic disks, and how it com- 
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1 Taking 12 + \og(0/H) = 8.66 for the solar metallicity (Asplund 
et al. 2005). 
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pares with more conventional galactic star formation has so far 
not been studied in detail. 

The number of reported "intergalactic" star forming re- 
gions has grown rapidly in the last few years, especially in 
groups (Mendes de Oliveira et al. 2004) and clusters of galaxies 
(Cortese et al. 2006). The area around NGC 5291 - an interact- 
ing system located at the edge of the Abell 3574 cluster - is 
one of the most spectacular objects where numerous instances 
of intergalactic HII regions are found. 

The system is composed of an early-type galaxy, 
NGC 5291, and a highly disturbed companion galaxy, "the 
Seashell". The first detailed observation of the system was car- 
ried out by Longmore et al. (1979) who detected on each side 
of NGC 5291 blue clumps identified as giant extragalactic HII 
regions. The spectra of some of them were obtained thanks 
to an object prism survey of HII galaxies (Maza et al. 1991; 
Pena et al. 1991). Malphrus et al. (1997) used the Very Large 
Array to map the neutral gas of the system. They discovered a 
huge asymmetrical ring-like structure connected to NGC 5291. 
It turned out to be one of the two most gas rich objects in the 
Southern hemisphere (Koribalski et al. 2004). The HI structure 
consists of a series of sub-condensations which correspond to 
the previously identified intergalactic HII regions. 

Due & Mirabel (1998) published a detailed optical (im- 
agery and spectra) and near-infrared study of the system. They 
found that the HII regions had moderately high metallicity 
(12 + log (O/H) = 8.4-8.6) and could thus exclude a pri- 
mordial origin for the HI structure. They claimed that the ring 
structure resulted from the past collision of an extended HI 
disk with a companion. They also showed that the most re- 
cent starbursts along the ring are younger than 5 x 10 6 years 
without any evidence for an underlying old stellar population. 
Braine et al. (2001) obtained with the SEST antenna CO(l - 0) 
and CO(2 - 1) spectra of two of the brightest condensations. 
Inferring the molecular gas content from the CO emission and 
comparing it with the SFR determined from Ha, they deter- 
mined that their star formation efficiency (SFE) was closer to 
that measured on average in spiral galaxies than in dwarfs of 
the same luminosity. Incidently, the detection of the millimetre 
lines was a further indication that the gas around NGC 5291 
had been previously enriched. Finally, Higdon et al. (2006) 
performed mid-infrared observations, both imaging and spec- 
troscopy, using the IRAC and IRS instruments on the Spitzer 
Space Telescope. They unambiguously detected Polycyclic 
Aromatic Hydrocarbon (PAH) features towards the two bright- 
est condensations of the system. Their spectral analysis also 
confirmed the young ages of the star forming episodes. 

In order to study in more detail the characteristics of the 
SF regions around NGC 5291, we collected ultraviolet images 
from the GALEX satellite, ground-based Ha images, and high 
resolution HI maps obtained with the Very Large Array which 
we combined with optical and mid-infrared Spitzer images of 
the system. See Figure 1 for a combination of those images. 

In section 2 we present the multiwavelength data set we 
have used. In particular, we describe how we defined and se- 
lected the individual intergalactic star forming regions and how 
we performed accurate aperture photometry on them. In sec- 
tion 3, we describe their main properties: overall distribution 



and star formation rates inferred from different tracers. We 
then compare them with those measured in various samples of 
galaxies and individual star forming regions in spiral disks. In 
section 4, we discuss the origin of the discrepancies between 
the various SF indicators and present the possible role of ex- 
tinction, metallicity and starburst age. We finally speculate on 
the possible impact of intergalactic star formation on its envi- 
ronment. 

Throughout this paper we usei/o = 72kms _1 Mpc~' which 
results in a distance for NGC 5291 of 62 Mpc. At this distance, 
l"=0.30kpc. 

2. Observations and data reduction 

2.1. Multi-wavelength data 

We describe in this section the various observations used for 
this paper. We present for the first time an ultraviolet map ob- 
tained with the GALEX space telescope. The high resolution 
HI maps were taken from the paper by Bournaud et al. (2007). 
The ESO/3.6m Ha data and Spitzer/IRAC mid-infrared im- 
ages were obtained from the literature or from the telescope 
archives. 

2.1 .1 . Ultraviolet images 

On May 2005 GALEX was used to acquire a set of NU V (near- 
ultraviolet; /i e ff = 227.1 nm and AA = 73.2 nm) broadband 
images of the target with a total exposure time of 2886 seconds. 
The field of view has a diameter of 1.24°, much larger than 
the NGC 5291 system. The PSF (point spread function) has 
a width of ~5.0". Due to technical difficulties at the time of 
observation, no far ultraviolet image (FUV) was obtained. 

2.1.2. Ha data 

Ha observations were carried out at the 3.6-m telescope of the 
European Southern Observatory (ESO) on April 2002, using 
the Fabry-Perot interferometer CIG ALE. The seeing was ~ I". 
The zeroth moment maps were extracted from the datacubes. 
For further technical details about the observations and data 
reduction, see Bournaud et al. (2004). 

Due to the small field of view of the camera, three individ- 
ual images were acquired, covering respectively the area cen- 
tred on the most luminous HII region to the south of NGC 529 1 , 
the central area close to NGC 5291 and the Seashell galaxies, 
and the area centred on the most luminous HII region to the 
north of NGC 5291. The southern and central images overlap 
partially but there is a gap between the central and northern 
images, as shown in Figure 1 . 

2.1.3. Mid infrared images 

The mid-infrared broadband images have been retrieved from 
the Spitzer public archives. The observations were carried out 
on February 2004 with the Infrared Array Camera (IRAC) on 
the Spitzer Space Telescope at 3.6 /im, 4.5 /im, 5.8 yum and 
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8.0 j-im. The resolution is ~3". The data have been prepro- 
cessed using the Spitzer Science Center (SSC) pipeline version 
SI 1.0.2. Due to the limited field of view (about 5'), the final 
image is a combination of two different pointings (each using 
a Reuleaux dither pattern), separated by 180" (see Figure 1). 
For further details about Spitzer data processing see the Spitzer 
Observing Manual 2 and the paper by Higdon et al. (2006). 

2.1 .4. HI observations 

The field containing NGC 5291 was observed with the NRAO 3 
Very Large Array (VLA) in the 21 -cm line of atomic neutral 
hydrogen. The observations and data reduction are fully de- 
scribed elsewhere (Bournaud et al. 2007). Here we will only 
dwell on the characteristics of the maps that were used in the 
analysis in this paper. 

We covered the entire velocity range of the system at 
48.8 kHz (~ 10.45 km s~'), after online Hanning smoothing. 
The AIPS 4 package was used to calibrate and Fourier transform 
the data. For this paper we used natural weighting which results 
in maps at 9.4 x 5.6 arcsec 2 resolution (and a beam position an- 
gle of 13°). The 1 <t rms noise is 0.33 mJy (corresponding to 
3.8 K). 

As the emission in this system is extended over more than 
lOarcmin, some flux due to a smooth, extended HI compo- 
nent might be missing from the interferometer data. It should 
be stressed that for spatial frequencies corresponding to scales 
smaller than ~ 5', all HI is recovered in our maps and that any 
contribution from an extended component becomes negligible 
in the small regions studied in this paper. 

2.2. Aperture photometry 
2.2.1 . Selection of the regions 

To identify the intergalactic star-forming regions belonging to 
the NGC 5291 system and reject from the subsequent analysis 
unrelated objects we used a combination of optical (B, V and 
R bands from Due & Mirabel 1998), Ha Fabry-Perot, mid in- 
frared (3.6, 4.5 and 8.0 pm) and HI observations. We used the 
following criteria in this order: 

- detection of Ha line emission at the NGC 5291 velocity 
on the first moment map of the Fabry-Perot datacube (see 
Bournaud et al. 2004), 

- proximity to the HI ring-like structure around NGC 5291, 
as this is the gas reservoir for star formation, 

- presence of diffuse emission in the near-ultraviolet map, 

- presence of extended, diffuse, blue regions on the true- 
colour BVR optical images, allowing the identification and 
subsequent removal of bright stars and distant background 
galaxies, 

2 http://ssc.spitzer.caltech.edu/documents/som/ 

3 The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 

4 The Astronomical Image Processing System is a software package 
provided by the National Radio Astronomy Observatory. 



- further rejection, based on their colours, of stars and back- 
ground galaxies using the image obtained combining the 
three mid-infrared bands. 

Because the maps obtained in the various wavelength bands 
do not cover completely the same area, we could in some re- 
gions only apply the HI and ultraviolet criteria. The Ha map, 
which of course provides an unambiguous detection of HII re- 
gions (apart from spurious lines falling right into the Ha red- 
shifted wavelength), unfortunately does not cover all of the HI 
structure. 

As HII regions are irregularly shaped and the system is 
crowded, we have chosen to define polygonal apertures. We 
could thus simultaneously take into account the variation of the 
PSF from one image to another avoiding as much as possible 
pollution from stars and from nearby objects, which would not 
have been easily possible using circular apertures. The shape 
and angular sizes of the polygonal apertures are the same for 
all bands and thus cover exactly the same physical region. 

Based on the above mentioned criteria, we have selected 
29 regions across the field, which correspond to intergalactic 
star-forming regions. They are indicated on Figure 3. 

For reference, we have also measured the fluxes towards 
the parent galaxies (NGC 5291 and the Seashell) and towards 
an HI peak just North of NGC 5291 with only faint ultraviolet 
and mid-infrared counterparts (information on Ha is lacking 
there). 

2.2.2. Measurements 

First of all, we have substracted the background in each band. 
To do so, we have manually measured the background level on 
each image using the IRAF 5 (Image Reduction and Analysis 
Facility) imexa procedure. Depending on the image, a single 
background was determined in case the background level was 
flat throughout the image (the ultraviolet image falls in this cat- 
egory as the standard pipeline already performs a background 
subtraction). Individually determined background levels had to 
be used for each of the regions in the Ha and infrared images. 
Each level is evaluated as the mean of several (up to 30) mea- 
surements of the background at different locations around the 
selected regions. The uncertainty on the background is the stan- 
dard deviation. 

We have also tested background removal using SExtractor 
(Bertin & Arnouts 1996). However, the large scale diffuse lu- 
minosity made the determination of the level quite uncertain 
and a small variation of the background extraction parameters 
caused a large variation of the flux, especially for the faintest 
regions where the difference between two extractions could be 
up to a factor of 3. 

The standard IRAC mid-infrared pipeline calibrates the im- 
ages in such a way that the flux of a point source is correctly 
determined when using a 10 pixel radius circular aperture and 
calculating the background in an annulus with inner and outer 

5 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 1. Composite image of the NGC 5291 system in pseudo-colours: 8.0 /im IRAC (red channel), Ha (green channel) and near ultraviolet 
GALEX (blue channel). The overplotted red and green contours represent respectively the 8.0 /im and the Ha fields of view. The blue contours 
represent the HI emission from the VLA HI map. The levels represent HI column densities of: 5.0 x 10 20 , 1.1 x 10 21 , 1.7 x 10 21 , 2.3 x 10 21 , 
2.9 x 10 21 ,3.5x 10 21 cm 2 . 



radii of 10 and 20 pixels 6 . This is due to the fact that the PSF 
is quite large, especially in the 8.0 /mi band where a significant 
part of the flux is scattered on arcminute size scales, which is 
much larger than 10 pixels. The calibrated units are set in such a 
way that the off aperture flux loss is exactly compensated when 
it is determined with a 10 pixel circular aperture. Therefore, 
when the aperture is larger (smaller) than 10 pixels of radius, 
the flux density is overestimated (underestimated). However in 



See http ://s sc . spi tzer. caltech. edu/irac/dh/ 



this paper we do not use a circular aperture for the determi- 
nation of the fluxes. Instead, the areas of interest are enclosed 
using polygons. As a result another method needs to be applied 
in order to estimate the flux outside our apertures and perform 
the necessary corrections. To do so we use the standard correc- 
tion values published in the IRAC manual for a point source 
correction, taking an equivalent radius of R= VS/tt, S being 
the surface in pixel units. The equivalent background annulus 
was taken with a radius from 10 to 20 pixels. Corrections for 
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Fig. 2. The different tracers of star formation towards NGC 5291: VLA HI map (top left), GALEX near-ultraviolet (top right), Fabry-Perot 
Ha (bottom left), and Spitzer 8.0 /im (bottom right) maps, all at the same spatial scale. The scaled 3.6 /im emission was subtracted from the 
8.0 nm map to remove contaminating stellar light from the mid-infrared emission (see for instance Pahre et al. 2004) and thus mainly shows 
dust emission. 



aperture areas falling between the published values have been 
interpolated. The uncertainty on the aperture correction there- 
fore depends on the aperture and the HII region shapes. The 
Poisson photon noise has also been taken into account. The 
absolute photometric calibration, according to the IRAC data 
handbook, is better than 5%. 

As the Fabry-Perot observations were not flux calibrated, 
we made an indirect calibration of the Ha data using the optical 
multi-object spectroscopic observations taken on July 1994 on 
the NTT ESO telescope at La Silla (Due & Mirabel 1998). The 
calibration is accurate to only about 30% due to uncertainties 
on the slit positions. The overlap of the two southern fields of 
view has been used to check the consistency of the calibration 
between the individual Fabry-Perot pointings. 



The absolute photometric uncertainty of the GALEX near 
ultraviolet image 7 is 7%. 

The fluxes derived in each band - when available - are 
presented in Table 1. The mid-infrared fluxes obtained have 
been compared whenever possible with the ones published by 
Higdon et al. (2006). The integrated fluxes of NGC 5291, the 
Seashell and PGC 048877 are consistent to within a few per- 
cent. Our photometry of individual star-forming regions agrees 
with the earlier values to within ~30% for the case of bright iso- 
lated regions, but it deviates substantially for the fainter ones. 
Our analysis indicates that the differences mainly originate 
from the way apertures were selected and how the background 

7 See Morrissey et al. (2005) and the GALEX Observers Guide 
http://www.galex.caltech.edu/DATA/grl_docs/ 
GRl_Observers_guide_vl.htm for more details. 
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Fig. 3. Ultraviolet image of NGC 5291 with a mosaic of the selected regions. The individual regions and the main galaxies in the group are 
labelled. The thumbnails are zoomed images of each selected region using pseudo-colours (red representing 8.0 /im emission, green Ha and 
blue ultraviolet), with the corresponding polygonal aperture overplotted in white. Note that for region 14 only the UV emission is available and 
shown (see Fig. 1). 
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substraction was performed: automatically by SExtractor in the 
Higdon et al. study, manually in this paper to make sure that re- 
gions with the same physical origin were enclosed. 

3. Results 

3.1. Distribution of the star forming regions 

The individual HI, NUV, Ha and 8.0 /im maps are presented in 
Figure 2 and a combination of them in Figure 1 . They show a 
remarkably similar morphology, indicating that all of them are 
physically linked and associated with the same process, most 
probably star-formation. We observe that star-forming knots 
are spread over ~180 kpc, all along the HI structure. They are 
actually located on or close to individual HI peaks. 

The brightest regions (in Ha, NUV and 8.0 /im) are located 
near the apparent tip and points of maximum column density 
of the HI ring: they correspond to our regions 1, 3 and 21; the 
latter two, also known as NGC 529 IS and NGC 529 IN, have 
been discussed in the literature as the possible progenitors of 
dwarf galaxies, and are probably associated with gravitation- 
ally bound gas clouds (Bournaud et al. 2004). As shown in 
Figure 2, the HI condensation just North of NGC 5291, which 
actually has one of the highest HI column density in the en- 
tire system (N(HI) > 3.3 X 10 21 cm~ 2 , hereafter referred to 
as HI concentration) shows only faint diffuse ultraviolet and 
mid-infrared counterparts 8 . The UV/HI and IR/HI flux ratios 
measured within this condensation deviate by an order of mag- 
nitude from the values determined from the other regions. One 
possibility would be that the starburst is extremely young and 
whatever O and B stars formed there have not dispersed their 
surrounding dense cloud material yet. Having a CO spectrum 
of that region and a far-infrared/millimetre image would help 
in checking this hypothesis. Having said that, inspection of in- 
dividual HI spectra across this cloud indicates that HI is found 
spread over a much wider range in velocities than at the other 
locations investigated. This could mean that what we are seeing 
in Figure 2 is a blend of clouds along the line of sight, adding 
up to result in an apparent massive HI concentration whereas 
each contributing cloud is much less massive, one of which be- 
ing just massive enough to have started to form stars. 

Figure 3 presents combined NUV, Ha and 8.0 //m emis- 
sion of the individual star-forming knots. The false-colour im- 
ages, all displayed with the same intensity scale, indicate that 
although they have the same morphology in each band, they 
show large variations in their relative fluxes. Local variations 
of the properties of the star forming regions are thus important, 
despite the fact that they all belong to a single structure proba- 
bly formed simultaneously during a collision. 

Another region, labelled number 1 in Figure 3, also shows a 
level of activity much higher than what would be expected from 
its HI column density. However the HI distribution of this ob- 
ject shows clear distortions: a gaseous tail without any 8.0 yum 
counterpart is seen just to the West of the star-forming region 

8 This region falls outside the Ha field of view. In the optical the 
region is also faint, the surface brightness averaged over the conden- 
sation being ~25.2 mag/" 2 in the B band. The global morphology of 
this region is similar to what is observed in the ultraviolet. 



(see Figure 1). There is evidence that in addition to the gravi- 
tational effects of the collision, the objects in the system may 
have suffered from the interaction of their interstellar medium 
with the intergalactic medium belonging to the cluster through 
which they move. A systematic offset between UV/Ha/MIR 
and HI peaks is observed with the HI being shifted to the oppo- 
site side of the cluster centre, which is located to the East of the 
galaxy. Why this process would be more efficient for object 1, 
which belongs to the same global structure as the other objects 
is however unclear. 

3.2. Star formation rates in the intergalactic medium 

The ultraviolet, Ha and infrared 8.0 fim emission are all known 
to be tracers of the level of star forming activity. The reliabil- 
ity of measuring the star formation rate based on each of these 
indicators, or a combination of them, has been extensively de- 
bated in the literature. It seems to depend very much on the type 
of environments as many external factors may affect it: contam- 
ination by old stars, level of dust extinction, etc. These tracers 
have not yet been calibrated for the special environments of 
collisional debris and IGM. We will thus, as a first step, use the 
classical calibrations obtained for spiral disks to get estimates 
of the SFR. From a comparison between the various estimates, 
we will then both discuss the validity of such tracers and see 
whether some global properties, such as the starburst age may 
be constrained. 

UV radiation is mainly emitted by massive stars younger 
than 100 x 10 6 years. The star formation rate has been 
calculated using the standard Kennicutt (1998) relation: 
SFR(UV)= [NUV] = 1.4xlO" 21 L v [WHr 1 ] M yr" 1 . This es- 
timator has been built under the assumptions of continuous star 
formation, with a Salpeter initial mass function (IMF; Salpeter 
1955) and mass cut-offs from 0.1 to 100 Mq. It is sensitive to 
extinction which reprocesses a fraction of the UV emitted by 
young stars which is reemitted in the infrared. 

Ionizing radiation (A < 91.2 nm) is mainly emitted by 
very massive, short-lived young stars. This radiation ion- 
izes surrounding gas clouds, which reemit partly in Ha, thus 
tracing starbursts. We have used again Kennicutt (1998): 
SFR(Ha)=[Z/a] = 7.9 x 10~ 35 L(Ha) [WJMoyr 1 , with the 
same IMF as used previously. This estimator is sensitive to ab- 
sorption of the ionizing photons by the molecular clouds sur- 
rounding young stars and thus very dependent on geometry 9 . 

At 8.0 //m we observe a combination of PAH broad emis- 
sion bands, continuum emission from hot VSG (very small 
grains) and stellar continuum. The photospheric stellar con- 
tribution at 8.0 fim is lower than 10% in our special environ- 
ment which is mostly devoid of an old stellar component (see 
below). The PAH are heated stochastically by far UV radia- 
tion from massive young stars and as such trace star forma- 
tion. Wu et al. (2005) built a SFR estimator based on a corre- 
lation between the 8.0 fim dust emission (stellar contribution 
being substracted using the 3.6 yum band) and radio or Ha lu- 

9 The radiation shoreward of 91.2 nm is more affected by dust ex- 
tinction than the near ultraviolet, therefore Ha is more sensitive than 
ultraviolet to the effects of geometry. 
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Table 1. Fluxes of selected regions in the ultraviolet, Ha, and mid-infrared (3.6, 4.5, and 8.0 /im). The "-" symbol indicates that the region 
falls outside the field of view, "<" means the aperture is heavily polluted by stars and the flux is therefore likely being overestimated and ">" 
means the aperture does not contain all of the object, e.g., to avoid a foreground object, the actual flux then being underestimated. 



Region 


GALEX name 


DM98 


F NUV 




F IR3.6 


F IR4.5 


F IR8.0 




(IAU name) 




G«Jy) 


(lO^WnT 2 ) 


0"Jy) 


0"Jy) 


0"Jy) 


1 


J134717.57-302821.ll 




274±20 


33.2+10.0 


331+101 


254+35 


640+99 


2 


J134722.91-302750.19 




92±7 


10.7+3.3 


144+41 


98+23 


629+41 


3 


J134723.00-302724.31 


a 


150+11 


18.8+5.7 


385+24 


251+32 


2088+55 


4 


J134719.76-302648.77 




53±4 


4.3+1.3 


131+12 


118+15.4 


207+27 


5 


J134723. 86-302646.25 


c 


98±7 


5.9+1.9 


221+28 


179+35 


805+63 


6 


J134720.69-302642.66 




14±2 


1.5+0.5 


31+5.4 


28+7 


97+12 


7 


J134724.51-302635.91 




14±1 


0.6+0.2 


53+12 


44+5 


73+12 


8 


J134717.79-302634.82 




6±1 


1.0+0.3 


31+7 


31+3 


60+7 


9 


J134723.95-302622.21 


d 


23±2 


1.5+0.5 


84+15 


61+6 


176+15 


10 


J134726.10-302602.63 




14±1 


1.4+0.5 


112+26 


95+10 


278+26 


11 


J134726.66-302552.06 


e 


34±3 


2.5+0.8 


199+18 


138+8 


384+18 


12 


J134726.45-302545.78 




31±3 


4.3+1.3 


229+34 


169+13 


509+34 


13 


J134723.37-302524.51 




12±1 


1.3+0.4 


188+18 


135+8 


311+19 


14 


J134700.56-302518.62 




62±5 


- 


- 


- 


- 


15 


J134727.77-302504.89 




5±1 


0.8+0.3 


70+9 


50+5 


189+10 


16 


J134717.92-302158.54 


f 


14±1 


1.4+0.4 


50+10 


41+6 


286+12 


17 


J134718.08-302150.30 




12±1 


1.3+0.4 


52+9 


37+6 


377+11 


18 


J134722.47-302144.13 


g 


30+2 


3.3+1.0 


146+28 


104+13 


1047+30 


19 


J134720.94-302131.97 




18±1 


1.6+0.5 


96+18 


51+7 


267+18 


20 


J134717.91-302121.48 


h 


87±6 


7.8+2.4 


287+58 


244+22 


1428+58 


21 


J134720.86-302055.70 


i 


259+18 


46.3+13.9 


638+38 


534+51 


5236+90 


22 


J134718. 18-302042.76 




30+2 


2.1+0.7 


63+27 


34+16 


427+27 


23 


J134721.66-302041.66 




17+1 


1.5+0.5 


14+17 


no detection 


114+17 


24 


J134717.07-302035.59 




21+2 


1.9+0.6 


1.7+20 


25+11 


31+20 


25 


J134717.97-302031.99 




15+1 


1.2+0.4 


38+13 


29+8 


235+14 


26 


J134721. 14-302026.07 




9+1 


0.7+0.2 








27 


J134719.74-302025.73 


j 


44+3 


4.2+1.3 


116+23 


100+14 


614+24 


28 


J134717.69-302018.20 




33+3 


2.3+0.8 








29 


J134717. 81-301958. 85 


k 


41+3 


2.5+0.9 








Seashell 


J134723. 17-302504.23 




98+7 


5.6+1.9 


14412+84 


8734+40 


5319+78 


NGC 5291 


J134717.81-302424.87 




>293±24 


>26.2±9.0 


>48619±385 


>29195±141 


>24846±374 


HI concentration 


J134721.02-302311.33 




71+8 




< 704+ 177 


<618±60 


1520+179 



minosities. The star formation rate calibration derived by them 
using the 8.0 yum/radio correlation is: SFR(IR)=[8.0] = 1.88 X 
10" 36 vL v (8.0/fln)[W]M Q yr- 1 . The Chary-Elbaz templates 
(Chary & Elbaz 2001) and the Kennicutt (1998) estimator 
based on the total infrared luminosity produce consistent re- 
sults. 

The star formation rates have been calculated using the re- 
lations given above and the fluxes presented in Table 1. We 
have corrected the fluxes for galactic extinction both in UV 
and in Ha using the Cardelli extinction curve (Cardelli et al. 
1989), and taking E(B-V)=0.064 from NED (NASA/IPAC 
Extragalactic Database). The systematic uncertainty of each 
estimator has not been taken into account. The UV SFR, 
[Ha] I [NUV] and [8.0] / [NUV] for each region is presented 
in Table 2. 

We stress that those estimators are likely not very adapted 
for our case (we expect a very young starburst, see Higdon 
et al. 2006). However they give more intuitive values to study 
the variations of the ratio of the SFRs in two bands, from one 
region to another. We use SFR(UV) as the reference since it 
is the most direct estimator (unprocessed star light) and more 



importantly, this is the only estimator which samples the star 
formation emission along the whole HI structure. 

The star formation rate derived from the NUV summed 
along the entire HI structure but excluding that of NGC 5291 
and the Seashell galaxies amounts to 1.37 ± 0.10 Mq yr _1 . The 
SFR of NGC 5291 and the Seashell is only 0.35 M G yr-\ or 
~20% of the SFR in the entire field of view. Thus, star forma- 
tion in this system is currently mostly observed in the inter- 
galactic medium. 

The total star formation rate in the regions in which we have 
measured fluxes in all three bands is: 

- SFR(UV)=[JVW] = 1.31 ±0.10 Moyr 1 

- SFR(Ha)=[#a] = 0.71 ± 0.22 M^yr 1 

- SFR(8.0)=[8.0] = 0.54 + 0.02 M yr 1 

The star formation rate derived from the ultraviolet is thus 
almost twice as large as that derived from Ha which itself is 
30% larger than that derived using the mid-infrared. 

To compare the different SFR estimators, we plot 
[Ha] I [NUV] (Figure 4) and [8.0] / [NUV] (Figure 5) versus 
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Table 2. Star formation rates deduced from the ultraviolet, Ha, and mid infrared, in solar mass per year. The symbol ":" means the estimation 
is uncertain. 



Region 


SFR(UV) 


[Ha] /[NUV] 


[8.0] / [NUV] 


1 


0.247+0.018 


0.57±0.15 


0.09+0.01 


2 


0.083±0.006 


0.55±0.15 


0.25+0.02 


3 


0.144±0.010 


0.56±0.15 


0.47+0.03 


4 


0.047±0.003 


0.39±0.10 


0.14+0.02 


5 


0.088±0.006 


0.29±0.08 


0.30+0.03 


6 


0.012±0.001 


0.53±0.17 


0.26+0.04 


7 


0.012±0.001 


0.21±0.06 


0.19+0.03 


8 


0.005±0.000 


0.86±0.24 


0.39+0.06 


9 


0.020±0.002 


0.31±0.09 


0.28+0.03 


10 


0.013±0.001 


0.46±0.15 


0.70+0.09 


11 


0.028±0.002 


0.34±0.10 


0.40+0.03 


12 


0.028±0.002 


0.65±0.18 


0.58+0.06 


13 


0.011 ±0.001 


0.51±0.14 


0.93+0.10 


14 


0.056±0.004 


- 


- 


15 


0.004±0.000 


0.79±0.27 


1.42+0.16 


16 


0.012±0.001 


0.48±0.12 


0.75+0.06 


17 


0.011±0.001 


0.50±0.14 


1.10+0.08 


18 


0.027±0.002 


0.52±0.14 


1.25+0.10 


19 


0.016±0.001 


0.43±0.12 


0.55+0.06 


20 


0.078±0.006 


0.43±0.12 


0.59+0.05 


21 


0.233±0.017 


0.85±0.22 


0.73+0.05 


22 


0.027±0.002 


0.33±0.10 


0.51+0.05 


23 


0.015±0.001 


0.41±0.12 


0.24+0.04 


24 


0.019±0.002 


0.42±0.12 


0.05+0.03 


25 


0.013±0.001 


0.38±0.11 


0.57+0.06 


26 


0.007±0.001 


0.39±0.10 




27 


0.040±0.003 


0.45±0.12 


0.40+0.04 


28 


0.030±0.002 


0.33±0.10 




29 


0.037±0.003 


0.29±0.09 




Seashell 


0.088±0.007 


0.27±0.11 


1.96+0.18 


NGC 5291 


>0.264±0.021 


0.42+0.18: 


3.05+0.29: 


HI concentration 


0.064±0.007 




0.77+0.17 




0.0 1 

0.001 



0.010 0.100 
SFR(UV)(M ffl .yr') 



0.0 1 

0.001 



h45 



h48 



47 

dl3 



si 



0.010 0.100 
SFR(UV) (Mo.yr 1 ) 



Fig. 4. Ratio of [Ha] / [NUV] versus SFR(UV). 



Fig. 5. Ratio of [8.0] / [NUV] versus SFR(UV). 



SFR(UV). On average, we get [Ha] / [NUV]=0A7 ± 0.16 and 
[8.0] / [NUV] = 0.53 ± 0.35. 

Figures 4 and 5 indicate strikingly large variations of 
[Ha] I [NUV] and [8.0] / [NUV] from one region to an other, 
even among the most luminous ones. For instance, there is a 



factor ~2 increase in [Ha] / [NUV] when comparing region 20 
with 21, or a factor ~6 decrease in [8.0] / [NUV] when going 
from region 1 to 3. Variations can also be important within a 
unique region, see for instance region 1 in Figure 3. As dis- 
cussed in Section 4, local variations of the star formation his- 
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tory, dust extinction and metallicity can in principle explain this 
scatter. This will be investigated further in the discussion. 

3.3. Intergalactic versus galactic star formation 

The intergalactic star-forming regions we have studied are ob- 
jects intermediate between isolated star-forming dwarf galaxies 
and individual HII regions in spiral disks. Indeed, the most lu- 
minous of them have simultaneously the structural properties 
and optical colours of typical primordial blue compact dwarf 
galaxies and the interstellar medium content, in particular the 
metallicity, of their parent disks. It therefore makes sense to 
compare the individual SF regions around NGC 5291 with both 
samples of individual galaxies, studied globally, and samples of 
individual HII regions within nearby galaxies (M51, M81 and 
Arp 82). Fortunately, several Spitzer based papers have recently 
provided the data required for these comparisons. 

We first compare the properties of the intergalactic star- 
forming regions in NGC 5291 with those measured globally 
for star-forming galaxies. For this comparison, we used three 
samples: the one studied by Rosenberg et al. (2006), consisting 
of 19 galaxies with a B-band absolute magnitude M B > -18.0; 
the sample of galaxies, mostly spirals, in the Spitzer SINGS 
legacy survey (Dale et al. 2007); and the Engelbracht et al. 
(2005) sample consisting of a heterogeneous collection of star 
forming galaxies without strong active nuclei including metal 
deficient blue compact dwarf galaxies. In Figure 6, we plot the 
8.0 pm fluxes, normalised to the 4.5 pm fluxes (a tracer of the 
hot dust emission) as a function of the monochromatic lumi- 
nosity at 8.0 pm. The 4.5 pm band was corrected for contam- 
ination by the photospheric stellar emission determined from 
the 3.6 pm band (which we assumed not to be significantly pol- 
luted by hot dust): Fi R 4 5 - aF lR i, 6, where a is a scaling coeffi- 
cient. We adopted a = 0.57 as used in Engelbracht et al. (2005). 
This is derived from Starburst99 models (Leitherer et al. 1999) 
which give typical values between 0.53 and 0.61; Engelbracht 
assigned a an uncertainty of 7%. 

Most of the individual regions of NGC 5291 exhibit 
roughly the same range of monochromatic luminosities at 
8.0 pm as the Rosenberg et al. sample of dwarf galaxies ex- 
cept for the brightest ones (see Figure 6). Based on their 
SFRs only, each of the 29 individual regions could by itself 
be considered a dwarf galaxy. Besides, they have on average 
a Fir8.o/(Fir4.5 - ffFiR3.6) ratio slightly, but not significantly, 
higher than the dwarf galaxy sample (see Table 3). The ratio for 
the SINGS is slightly larger but very similar. The Engelbracht 
et al. sample, in spite of its inhomogeneity, shows a similar 
mean value as for the NGC 5291 regions. 

In figure 7 we compare [NUV] to [8.0] for the galax- 
ies of the SINGS sample and the HII regions detected in the 
NGC 5291 system. 

We notice that the least luminous galaxies from the SINGS 
sample present an important ultraviolet excess. We think it is 
due to the PAH defiency seen in the low metallicity galaxies, 
rather than an ultraviolet excess per se. On the other hand, the 
most luminous SINGS galaxies present a strong infrared ex- 
cess, due to the higher metallicity and the large population of 




°0 

□ o o 



□# o ° 



8i> 



10" 10" 
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Fig. 6. Ratio of the mid-infrared fluxes Fi R8 o/(Fir4 5 - aFiR 3 6 ) versus 
the monochromatic luminosity at 8.0 fim of SINGS galaxies (Dale 
et al. 2007, open circles), of the dwarf galaxies in the Rosenberg et al. 
sample (triangles), of Engelbracht et al. sample (squares) and the HII 
regions in NGC 5291 (filled circles with error bars). The large filled 
black circle represents the sum over all NGC 5291 regions. 
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Fig. 7. Star formation rate estimated from the IRAC 8.0 /im band as 
a function of SFR(UV) for the NGC 5291 intergalactic SF regions 
(filled circles) and for the galaxies of the SINGS sample (Dale et al. 
2007, open circles). 



non-ionizing stars which can excite the PAH in addition to mas- 
sive young stars. 

The relative scatter (see Table 3) for all these samples are 
quite similar. This is surprising as the Rosenberg et al., the 
Englebracht et al. and the SINGS samples are composed of het- 
erogeneous dwarf galaxies whereas the NGC 5291 HII regions 
lie in the same HI structure. 
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Fig. 8. Ratio of the mid-infrared fluxes Fir8.o/(Fir4.s - aFiR3.6) versus 
the monochromatic luminosity at 8.0 jim of M81 (Perez-Gonzalez 
et al. 2006, purple diamonds), Arp 82 (Hancock et al. 2007, red 
squares; filled squares represent star forming regions in tidal features) 
and NGC 5291 (filled circles with error bars) HII regions. The large 
filled circled represents the integration of all NGC 5291 regions. 

In Figure 8 we perform a similar comparison as in Figure 6 
but with HII regions in the nearby spiral M8 1 10 (for more de- 
tails see Perez-Gonzalez et al. 2006) and the interacting system 
Arp 82 (Hancock et al. 2007). 

The Firs.o/(Fir4.5 - aFny.e) ratio of the regions in M81 is 
roughly similar to that of the most luminous HII regions of 
NGC 5291 even though their 8.0 fim monochromatic luminos- 
ity is lower by one to two orders of magnitude. Beside, for 
the less luminous regions, the ratio is higher in M81 which 
could indicate a deficiency in PAH emission in the HII re- 
gions of NGC 5291. The star forming regions in Arp 82, es- 
pecially along the tidal arms of the interacting system, follow 
more closely what can be seen in NGC 5291 both in the 8.0 pm 
flux and the F IR8 /(F 1R 4 5 - aF IR 3 6). Therefore both systems 
have on average the same PAH properties compared to the dust 
continuum. 

In Figures 9 and 10, respectively, we present a comparison 
of the [NUV] to [8.0] and [NUV] to [Ha] SFRs for individual 
HII regions in the HI intergalactic structure of NGC 5291, for 
the HII regions in the disk of M51 (Calzetti et al. 2005), M81 
(Perez-Gonzalez et al. 2006) and Arp 82 (Hancock et al. 2007). 

We noticed in section 3.2 that the SFRs of the NGC 5291 
HII regions estimated from the UV fluxes were on average 
higher by a factor of 2 than those obtained from the Ha and 
8.0 fim bands 1 1 . Figures 9 and 10 show that this systematic UV- 

10 The mid-infrared flux densities were provided by Perez-Gonzalez 
(priv. comm.). 

11 Using the standard calibration of Kennicutt (1998) and the as- 
sumptions made therein. 



Fig. 9. Star formation rate estimated from the IRAC 8.0 fim band as a 
function of SFR(UV) for the NGC 529 1 intergalactic SF regions (filled 
circles), for samples of HII regions in two spiral disks: M5 1 (Calzetti 
et al. 2005, blue triangles; filled triangles represent inner regions 
where Paa flux is available) and M81 (Perez-Gonzalez et al. 2006, 
purple diamonds), and for the interacting system Arp 82 (Hancock 
et al. 2007, red squares; filled squares represent star forming regions 
in tidal features). 

excess is not observed for the individual HII regions of M5 1 
and M81; the regions in M51 present an IR excess, whereas 
those of M8 1 seem to be too Ha rich. No systematic UV ex- 
cess can be detected in the Arp 82 interacting system. 

More quantitatively, when normalised to the 8.0 fim band, 
the NUV flux in NGC 5291 is in excess by a factor of 8, 2 and 
1 .6 compared to M5 1 , M8 1 and Arp 82 respectively. When nor- 
malised to Ha, the corresponding excess factors are 2, 6 and 2. 
Finally, one should note that NGC 5291 presents a scatter in 
the Ha vs. UV relation which is smaller than for the compari- 
son galaxies, while the MIR vs. UV is roughly comparable. 

In summary, while globally the 8.0 fim emission in the in- 
tergalactic star forming regions of NGC 5291 are not signifi- 
cantly different than those of star-forming dwarf galaxies, they 
seem to present with respect to individual HII regions in spirals, 
and with the Ha and MIR SF indicators, a clear UV excess. 

4. Discussion 

In the previous section, we presented the properties of the star 
forming regions along the huge HI structure around NGC 5291, 
a study based on three SF indicators: UV, Ha and the mid- 
infrared. We found that the three of them were correlated, as for 
more classical galactic star-forming regions, but showed with 
the latter some differences, in particular a UV excess. We dis- 
cuss here the origins of these differences and of the scatter in 
the correlations. The goal is to determine which SF indicator 
- or combination of them - may be used to estimate the total 
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Table 3. Mean flux and SFR ratios for the different samples. Extremely discrepant points (if any) have been discarded before performing the 
calculation. 



Sample 


(FlR8.o/(FlR4.5 - QrF 1R3 6 )> 


([S.0]/[NUV]) 


([Ha] 1 [NUVfi 


NGC 5291 


21 ± 17 


0.5 ± 0.4 


0.5 ±0.2 


Rosenberg 


17 ± 12 






SINGS 


25 ± 16 


3.3 ±7.0 




Engelbracht 


24 ±23 






M51 




4.2 ±5.5 


1.1 ± 1.3 


M81 


50 ±35 


1.2 ±0.9 


2.8 ±2.2 


Arp82 


45 ±24 


1.8 ± 1.8 


1.0 ±0.7 



1.0000 



0.1000 r 



g 0.0100- 



a 



0.0010 



0.0001 r 




0.0001 



0.0010 0.0100 

SFR(NUV)(M s .yf') 



0.1000 



1.0000 



Fig. 10. Star formation rate estimated from the [Ha] luminosity as 
a function of SFR(UV) for the NGC 5291 intergalactic SF regions 
(filled circles), for samples of galactic HII regions in two spiral disks: 
M51 (Calzetti et al. 2005) (blue triangles; filled triangles represent 
HII regions near centre of M5 1 where Pao- flux is available) and M8 1 
(Perez-Gonzalez et al. 2006, purple diamonds), and for the interact- 
ing system Arp 82 (Hancock et al. 2007, red squares; filled squares 
represent star forming regions in tidal features). 



SFR in this special environment, bearing in mind that the previ- 
ous calibrations were derived for galactic environments. Once 
a "total" SFR is determined, we briefly discuss its impact on 
the IGM. 



4.1. The effect of metallicity and PAH strength on the 
MIR star formation indicator 

In most galactic star-forming regions, the dust emission around 
8.0 fim is largely dominated by PAH features. The latter 
are however known to be absent in metal-poor objects such 
as the Blue Compact Dwarf Galaxies (Galliano et al. 2005; 
Engelbracht et al. 2005; Madden et al. 2006; Wu et al. 2006). 
Thus for these systems, the IRAC based mid-IR-SFR cali- 
bration, as the one used in this study, could be fairly inac- 



curate. The intergalactic HII regions around NGC 5291, al- 
though they share many properties with those of BCDGs, lie 
however within an HI structure which had been previously pre- 
enriched. Due & Mirabel (1998) measured an oxygen abun- 
dance of 12 + log(0/H) = 8.4 to 8.6, with an average of 
8. 49^'^. Although this is somewhat below the solar value, 
this is consistent with a PAH detection. Indeed, the IRS spec- 
tra obtained by Higdon et al. (2006) indicate that the IRAC 
8.0 fim emission of the two most luminous star forming regions 
is mainly due to PAHs. 

In Figure 1 1, we plot the PAH emission normalised to the 
dust continuum as a function of metallicity for the NGC 5291 
star forming regions, the dwarf galaxies from the Rosenberg 
sample and the galaxies from the Engelbracht sample. The lu- 
minosity at 8.0 fim appears to be consistent with that of dwarf 
galaxies with the same metallicity. Wu et al. (2006) directly 
measured a correlation between the equivalent widths of the 
PAH features and metallicity. We checked this for region 21. 
We found that the EWs measured in the IRS spectrum pub- 
lished in Higdon et al. (2006) 12 , using a method similar as in 
the Wu et al. (2006) paper, were consistent with the Wu et al. 
correlation. 

Thus, on average, the relative strength of the 8.0 fim emis- 
sion in NGC 5291 appears to be normal, i.e. consistent with 
that expected for star-forming regions with the same metallic- 
ity. Therefore the calibration derived between the IRAC mid- 
infrared emission and the SFR so far only validated for spiral 
disks should also apply to the collisional debris. 

The large scatter in the Firs ,o/(Fir4.5 - aF^^) flux ratio 
along the system (see Figures 6 and 8) is however surpris- 
ing. In any case, it cannot be simply explained by local vari- 
ations of the metallicity as the latter was measured by Due 
& Mirabel (1998) to be quite uniform all along the HI struc- 
ture. For instance, very bright regions, such as 1 and 3, show 
discrepant F^g o/(Fir4.s - aFiR3.6) flux ratios although they lie 
very close on the sky and thus should have very similar metal- 
licities. Unfortunately no optical spectrum of region 1 is avail- 
able to check this. The IRS spectra of the two most luminous 
regions (regions 3 and 21 in this paper) also show variations 
in their PAH strengths and dust continuum shortwards of 5 
fim and above 9 fim, which are not yet understood. Clearly, 
the metallicity is not the only parameter playing a role in the 



12 EW(6.2 yum)=0.58 fim, EW(7.7 yum)=0.52 fim, EW(8.6 fim)=0.53 
fim, EW(11.2 y um)=0.65 fim. 
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Fig. 11. Ratio of the mid-infrared fluxes Fi R8 .o/(Fir4.5 - aF\ Ri 6 ) (with 
a = 0.57) versus the metallicity of dwarf galaxies from the Rosenberg 
et al. sample (triangles), Engelbracht galaxies (open squares) and HII 
regions in NGC 5291 (filled circles). 



Fig. 12. [Ha] I [NUV] as a function of [8.0] / [NUV] for HII knots 
with an uncorrected ultraviolet star formation rate greater than 0.02 
Moyr~'. Filled circles are uncorrected for extinction while empty 
ones are. 



8.0 yum emission, as also shown in several studies of galactic 
star formation. 

4.2. The effect of dust extinction on the UV and Ha 
star formation indicator 

We have seen earlier that dust is undoubtedly present all along 
the HI ring. Could local variations of the dust extinction and ge- 
ometrical factors account for the large scatter in the individual 
SF properties? To examine that, we have corrected both [NUV] 
and [Ha] for extinction, the value of which was determined 
locally towards each condensation. The extinction is actually 
composed of two components: one due to the dust clouds sur- 
rounding newborn stars, which is due to the molecular gas and 
very dependent on geometry and may vary dramatically from 
one region to the other; a second one is due to the dust along 
the line of sight, which can be approximated by a dust screen. 
Assuming the ratio of dust to HI is constant, we can derive the 
extinction variation from the HI maps. 

To quantify this HI based extinction, we used previous ex- 
tinction estimates determined by Due & Mirabel (1998) from 
the measurement of the Balmer decrement in several associ- 
ated HII regions. Correlating the peak HI column densities and 
the optical extinction in the V band with 7 data points (after 
having excluded 2 which showed large discrepancies), we ob- 
tained a calibration of the Hi-based extinction. We estimate the 
uncertainty to be about 0.4 dex with this method. 

To convert the extinction from the V band to other bands, 
we used the average Large Magellanic Cloud A (A) /A v of 
Gordon et al. (2003). We find A NUV /A V = 2.80 and A Ha /A v = 
0.89. 

We have corrected the NUV and Ha fluxes for dust extinc- 
tion based on the calibrated HI column densities. Results of this 
correction are shown in Figure 12 for objects with uncorrected 
ultraviolet star formation rate greater than 0.02 Mq yr _1 , which 
have the most secure determination of the SFR. 



The relative scatter of [Ha] / [NUV] and [8.0] / [NUV] is 
not reduced by the extinction correction. Clearly, variations of 
the dust extinction alone cannot explain all the scatter in the 
SF properties found along the HI structure, unless very local 
variations dominate. 

4.3. The effect of starburst age on the UV and Ha star 
formation indicators 

The so-called "UV-excess" seems to be a determining char- 
acteristic of the intergalactic star forming regions around 
NGC 5291. This was also observed in a few other intergalactic 
SF regions (Werk et al. 2006). We discuss here whether this ex- 
cess, as well as the large scatter of the [Ha] / [NUV] after the 
correction for dust extinction, could result from some starburst 
age effects. 

Due & Mirabel (1998) measured H/3 equivalent widths as 
large as 140 A of several HII regions around NGC 5291, which 
is indicative of young starbursts (see for instance Terlevich 
et al. 2004). The detection of Wolf-Rayet features in the most 
luminous one shows that the current star formation episode 
there is less than 5 X 10 6 years old. More recently, Higdon et al. 
(2006) determined ages of ~ 5 x 10 6 years from [Ne II]/[Ne III] 
emission line ratios for the objects for which mid-infrared IRS 
spectra were available. Are these indications for starbursts be- 
ing young consistent with the UV excess observed in all in- 
tergalactic HII regions around NGC 5291? At first sight, one 
may think that a UV excess, or conversely, an Ha deficit would 
rather suggest older, declining bursts since the UV time scale 
is a factor of ten larger than the Ha one (see section 3.2). 

To pursue the analysis, we have modelled the HII regions 
using the spectral evolution code Pegase II (Fioc & Rocca- 
Volmerange 1997, 1999). We have used a set of parameters 
differing by the star formation history (SFH), going from an 
instantaneous starburst to exponentially decreasing ones with a 
timescale between 1 and 40 million years. We used a Salpeter 
IMF from 0.1 to 100 M and a constant metallicity Z=0.008 
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(65% solar metallicity, corresponding to 12 + log (O/H) = 
8.47), in agreement with the metallicity determined by Due & 
Mirabel (1998). Extinction from the gas has not been taken into 
account in the simulations in order to isolate age effects only. 
For comparison we therefore need to use the extinction cor- 
rected ratios. We have in particular studied how the Hfi equiv- 
alent width and the [Ha] / [NUV] ratio evolve with SFH (see 
Figures 13 and 14). The first plot suggests that the ongoing 
star formation may not have lasted for more than 10 x 10 6 
years for the regions presently showing the largest H/3 equiv- 
alent width, whatever the timescale assumed for the starburst. 
Other than that, very young, quasi-instantaneous starbursts ig- 
nited 3-5 xlO 6 years ago may well cohabit with older and more 
extended star-formation episodes (Neff et al. 2005), i.e. in re- 
gions showing smaller H/? equivalent width. The second plot 
shows that the observed range for [Ha] / [NUV] may be ob- 
tained for various SFH, but that in any case, all bursts are well 
advanced. A constant or just slowly declining SFR is incon- 
sistent with the data. However, our dataset may accommodate 
quasi-instantaneous bursts which occurred 4 - 10 x 10 6 years 
ago up to bursts with slightly longer time scales of 10 x 10 6 
years which started about 40 x 10 6 years ago. Combining the 
two criteria, we conclude that the regions that are currently the 
most active ones (i.e. region 1, 3 and 21) are all very young but 
should have had their SF peak already 5 x 10 6 years ago while 
the more quiescent ones may be a few 10 x 10 6 years older, still 
sustaining star formation but at a lower rate than before. 

One should note that quasi-instantaneous young starbursts 
for all HII regions around NGC 5291 are also fully compat- 
ible with the data, i.e. the H/3 equivalent width, the UV ex- 
cess and the somehow large region-to-region variations of the 
[Ha] I [NUV] to UV ratio. The latter may be understood noting 
how this ratio changes quickly after the beginning of the star- 
burst: a small age difference can produce a considerable differ- 
ence in [Ha] / [NUV] (see Figure 14). Probing the presence, or 
actually the absence, of an intermediate age stellar population 
of a few 10 x 10 6 years would be required to validate the sce- 
nario of a rapidly decreasing starburst. This requires a model of 
the full Spectral Energy Distribution of each region, from the 
UV to IR, which is beyond the scope of this study. This task 
is somehow easier in the case of the NGC 5291 intergalactic 
star-forming regions. Indeed, they should not contain a stellar 
component of a few 100 x 10 6 years or older coming from the 
parent galaxies. 

As all the HII regions seem to have roughly the same age 
it is legitimate to wonder how star-formation may have started 
quasi simultaneously along a structure as extended as 180 kpc, 
and how likely we are to observe it so close to its onset. 
One should however remember that the HI ring which hosts 
all these regions was probably formed during a single violent 
event, likely the collision of a previously existing extended HI 
disk with a companion galaxy (Due & Mirabel 1998). Given 
its ring-like shape and overall velocity field, the HI structure 
was most probably not formed by tidal forces but rather by 
an expanding density wave following a bull's eye collision. 
According to numerical simulations (Bournaud et al. 2007), 
the ring around NGC 5291 has an age between 300 x 10 6 
and 400 x 10 6 years, very different from the star formation 




Age [Myr] 

Fig. 13. Plot of the variation of the equivalent width of H/J in A, as a 
function of ages for different types of starbursts. The black line shows 
an instantaneous starburst whereas the red one corresponds to an ex- 
ponentially decreasing starburst with a timescale t = 40 x 10 6 years. 
Lines for all intermediate t with intervals of 1 x 10 6 years are also 
shown. The light grey zone indicates the range of observed values of 
the H/3 equivalent width for the few regions where this was measured 
(Due & Mirabel 1998). 




Age [Myr] 

Fig. 14. Plot of the variation of the [Ha] / [NUV] ratio as a function 
of age for different types of starbursts. The black line shows an instan- 
taneous starburst whereas the red one corresponds to an exponentially 
decreasing starburst with a timescale r = 40 x 10 6 years. Lines for all 
intermediate t with intervals of 1 x 10 6 years are also shown. The light 
grey zone indicates the range of observed values of [Ha] / [NUV] af- 
ter correction for extinction. 



timescale. However, in the case of a collisional ring and un- 
like tidal tails, the star formation can be suppressed during the 
first few 100 x 10 6 years, i.e., until the HI ring expansion stalls 
after which star formation proceeds. If that scenario is indeed 
correct, this would bring the dynamical and the star-formation 
time scales together and explain the, at first sight, rather intrigu- 
ing observation of quasi simultaneous star formation along the 
entire ring. 
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4.4. The total intergalactic star formation rate In 
NGC 5291: impact on the IGM 

The integrated star formation rate along the HI structure - un- 
corrected for dust extinction - derived from the GALEX NUV 
emission, amounts to 1.31+0.10 M yr" 1 , 0.71+0.22 M yr _1 
in Ha and 0.54+0.02 Mq yr" 1 from the 8.0 /itn band emission 
(see section 3.2). 

Whether the uncorrected UV, Ha and IR SFRs should be 
added to get the total, dust free SFR, or whether a simple ex- 
tinction correction for each tracer is enough is a matter of active 
debate (Hirashita et al. 2003). Because our estimates of the ex- 
tinction, based on the HI column density, are rather uncertain, 
we decided not to apply a dust correction for each individual re- 
gion. Instead we decided to add the uncorrected SFR obtained 
from Her to the SFR based on the dust emission scaled by a 
coefficient which takes into account the amount of ionizing ra- 
diation (actually radiation emitted by the same stars as traced 
through Hot) absorbed by PAH. The latter was determined em- 
pirically by fitting SFR(Ha)+kxSFR(IR) to the extinction cor- 
rected SFR(Ha) using only datapoints for which measures of 
Ay obtained from the Balmer decrement were available. The 
advantage of using a combination of Ha and IR rather than UV 
and IR (as commonly done) is that the former two indicators 
probe SF with relatively similar time scales of 10 x 10 6 years 
(rather than 100 x 10 6 years for UV). Doing so, we obtained 
k = 0.27 and a total SFR all along the HI structure of 0.83 
Moyr" 1 . This value is a factor of 1.7 smaller than the value 
determined from the uncorrected near UV emission, confirm- 
ing again the earlier mentioned UV excess. 

However, we should note at this stage that the SFRs were 
estimated using calibrations valid for a constant star formation 
over 100 x 10 6 years (Kennicutt 1998). As mentioned before, 
the star formation history in NGC 5291 may have been quite 
different. There are various independent indications pointing 
towards young, quasi-instantaneous, starbursts and the measure 
of the actual SFR may be different. 

For young and rapidly decreasing starbursts, models (see 
Figures 15 and 16) indicate that the Kennicutt UV and Ha 
calibrations overestimate by factors of 2 to 4 the actual SFR. 
Thus if the quasi-instantaneous starburst hypothesis is correct, 
the real "intergalactic" star formation rate around NGC 529 1 
should be 0.3-0.7 Mq yr" 1 , still higher than for typical dwarf 
and irregular galaxies (Hunter & Elmegreen 2004). However, 
one should keep in mind that the UV-excess indicates that the 
star formation rate was higher a few 10 6 years ago, with values 
close to that observed in quiescent spiral galaxies. 

We thus find that around NGC 5291, star-formation pro- 
ceeds at a rate typical of spirals but in an environment which 
is much less dense, and where the effect of stellar winds on 
the surrounding medium should be much more dramatic. As 
shown, among others, by Ferrara & Tolstoy (2000), the ability 
of an object to lose its heavy elements, and thus enrich the inter- 
galactic medium, depends very much on the mass and shape of 
the halo it sits in. Collisional debris being dark matter deficient, 
they should have a reduced capability to retain their metals. 
How the particularly active star-forming regions in NGC 529 1 
manage to survive this feedback of mechanical energy still re- 




Fig. 15. Plot of the variation of [Ha] (it is to say the star formation rate 
that would give the Kennicutt estimator) over the actual star formation 
rate as a function of age for different types of starbursts. The black line 
shows an instantaneous starburst whereas the red one corresponds to 
an exponentially decreasing starburst with a timescale t = 40 x 10 6 
years. Lines for all intermediate t with intervals of 1 x 10 6 years are 
also shown. 
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Fig. 16. Plot of the variation of the [NUV] (it is to say the star forma- 
tion rate that would give the Kennicutt estimator) over the actual star 
formation rate as a function of age for different types of starbursts. 
The black line shows an instantaneous starburst whereas the red one 
corresponds to an exponentially decreasing starburst with a timescale 
t = 40x 10 6 years. Lines for all intermediate r with intervals of 1 x 10 6 
years are also shown. 



mains to be studied. Given their large gas reservoir, it is some- 
what surprising that they are currently already fading, unless 
SF has ceased as a result of feedback from massive stars. 

Clearly, a system like NGC 5291 is exceptional in the lo- 
cal Universe. Studying a sample of interacting systems (Struck 
2006; Smith et al. 2007; Hancock et al. 2007) found that their 
external plumes and tails only account for at most 10% of the 
total star formation. In NGC 5291, from the UV band we can 
see that more than 80% of the stars are currently formed along 
the intergalactic ring. When using the 8.0 fim band, removing 
the old stellar content from the NGC 5291 galaxy (this popu- 
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lation is responsible for roughly 50% of the 8.0 pm flux), the 
intergalactic SF still accounts for more than 50% of the SFR. 
In the distant Universe, however, such large gaseous structures 
may have been more common 13 and the star formation regions 
they hosted may have contributed significantly to the direct en- 
richment of the IGM. 

5. Conclusion 

We have performed a multi-wavelength study of the NGC 5291 
HII regions using ground-based, GALEX, and Spitzer Space 
Telescope data of 29 intergalactic HII regions around 
NGC 529 1 . These regions lie along a previously identified huge 
HI ring-like structure of collisional origin. We have measured 
their fluxes in the near ultraviolet, Ha, and mid-infrared in or- 
der to compare these different star formation tracers. The mor- 
phology observed in the three bands is remarkably similar, in- 
dicating that each can serve as a star formation tracer. We have 
shown that the intergalactic IR emission at 8.0 pm, dominated 
by PAH, when normalised to the hot dust continuum, is compa- 
rable to the integrated emission of dwarf galaxies of the same 
metallicity and to the emission of individual HII regions in spi- 
rals. This indicates that the 8.0 pm band emission is probably 
an estimator of the SFR that is as reliable for these regions as 
it is for spirals. Whereas the 8.0 pm is "normal", there is a 
clear excess of near ultraviolet emission compared to individ- 
ual HII regions in spirals. The [8.0] / [NUV] and [Ha] / [NUV] 
are low although there are some strong variations from one 
region to an other. The variations of the metallicity sensitive 
[8.0] / [NUV] ratio cannot be explained by a metallicity effect 
since the oxygen abundance was found to be uniform along 
the entire HI structure. Moreover, correcting for the spatial 
variations of the dust extinction does not reduce the scatter in 
the [Ha] / [NUV] and [8.0] / [NUV] SFR ratios. Their varia- 
tions are best explained by age effects. A model of the evolu- 
tion [Ha] I [NUV] with time favours young instantaneous star- 
bursts: the ultraviolet excess or deficit of Ha plus the large 
scatter indicating rapid variations of the [Ha] / [NUV] ratio 
are difficult to explain with a constant SFR or moderately de- 
creasing SFR. A young starburst is consistent with constraints 
from the H/3 equivalent widths and mid-infrared Nell fluxes. So 
far no indication for the presence of an old stellar component 
has been found. A detailed study of the full SED of the star- 
forming regions, including dust models, would be required to 
confirm it. Comparing the individual starburst ages with the 
dynamical timescale for the formation of the HI ring in which 
they lie, we concluded that the apparent observation of a quasi- 
simultaneous onset of star formation over the 180 kpc long 
gaseous structure can be understood if SF commences only af- 
ter expansion of the ring has stalled, after several 100 x 10 6 yr. 

Using the standard Kennicutt calibrations, we estimated 
from the UV, uncorrected for dust extinction, a total SFR 
around NGC 5291 of 1.3 M©yr _1 . If the quasi-instantaneous 
starburst hypothesis is correct, then this current value of the 

13 The gigantic Lyo- blobs found in the distant Universe (Nilsson 
et al. 2006; Jimenez & Haiman 2006) actually present some similari- 
ties with NGC 5291. 



SFR may have been overestimated by factors of 2 to 4. In 
that case, the SFR would have been much higher in the recent 
past. Finally it appears that about 80% of star formation in the 
NGC 5291 system takes place in the intergalactic HII regions, 
a value much higher than in nearby interacting systems. 
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